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DOUBLE QUASARS: PROBES OF BLACK HOLE SCALING RELATIONSHIPS AND MERGER SCENARIOS 
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ABSTRACT 

We analyze the available sample of double quasars, and investigate their physical properties. Our 
sample comprises 85 pairs, selected from the Sloan Digital Sky Survey (SDSS). We derive physical 
parameters for the engine and the host, and model the dynamical evolution of the pair. First, we 
compare different scaling relationships between massive black holes and their hosts (bulge mass, 
velocity dispersion, and their possible redshift dependences), and discuss their consistency. We then 
compute dynamical friction timescales for the double quasar systems to investigate their frequency and 
their agreement with the "merger driven" scenario for quasar triggering. In optical surveys, such as the 
SDSS, N double, qso/N qso as 0.1%. Comparing typical merging timescales to expected quasar lifetimes, 
the fraction of double quasars should be roughly a factor of 10 larger than observed. Additionally, 
we find that, depending on the correlations between black holes and their hosts, the occurrence of 
double quasars could be redshift-dependent. Comparison of our models to the SDSS quasar catalog 
suggests that double quasars should be more common at high redshift. We compare the typical 
separations at which double quasars are observed to the predictions of merger simulations. We find 
that the distribution of physical separations peaks at ~ 30 kpc, with a tail at larger separations 
(~ 100 — 200 kpc). The peak of the distribution is roughly consistent with the first episode of quasar 
activity found in equal mass mergers simulations. The tail of the quasar pairs distribution at large 
separations is instead inconsistent with any quasar activity predicted by published simulations. These 
large separation pairs are instead consistent with unequal mass mergers where gas is dynamically 
perturbed during the first pericentric passage, but the gas reaches the black hole only at the next 
apocenter, where the pair is observed. 

Subject headings: quasars: general, galaxies: active, black hole physics, galaxies: evolution 
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1. INTRODUCTION 

Interactions and mergers are an integral part of the 
process that builds up the galaxy population. Mergers 
are expected to be the dominant growth route of galax- 
ies in hierarchic al scenarios linked t o Cold Dark Mat- 
ter cosmologies (jWhite k Reesl fl978l ) . Recent high res- 
olution observations have shown evidence of such phe- 
nomena, revealing that many of the most massive galax- 
ies at z ~ 1 — 4 are inde ed interacting/merging galax- 
ies (|Chapman et al.l 120031 e.g.,). Many of submillimc- 
ter detected galaxies also show signs of nuclear activ- 
ity (|Alexander et all I200H 12001 [20051) . and are found 
near radio galaxies, suggesting that quasar activity was 
common in dense environments at high redshift. The- 
oreti cal models indeed link quasar activity to mergers 
Ce.g.JHaehnelt fc Reeslll993l:ICavaliere k VittorinilpMl 
Kauffmann k HaehnelH 120001: IWvithe k Loebl \20o£ 
Volonteri et all 120031 : IDi Matteo et all I2005D . If most 



galaxies host a central massive black hole (MBH) , as im- 
plied by the local observations, and nuclear activity is 
triggered in mergers, double/binary quasars should also 
be common in a hierarchically evolving Universe. The 
expected fraction of double quasars expected in a hier- 
archical scenario w here mergers trigger nu clear activity 
was investigated by I Volonteri et all(|2003[ ). If accretion 
were assumed to be triggered on both MBHs involved in 
a major merger (that is, mergers between galaxies with 
similar masses, up to a factor of 10 difference), the frac- 
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tion of binary quasars would be about 10% , in conflict 
with the existing data (Vo lonteri et al. 2003) . A rough 
agreement with the data could be reached assuming that 
only the MBH hosted in the larger halo accretes and ra- 
diates during a major merger. In this scenario at least 
two subsequent major mergers are required to give ori- 
gin to a double quasar system, the satellite MBH being 
activated in the first merger event. In these models, the 
mean timescale over which MBHs accrete and radiate is 
of order ~ 4 x 10 7 years, much short er than the typ- 
ical t ime between two major mergers (|Somerville et all 
12001 . 

Only 16 pairs are known in the Large Bright Quasar 
Survey ijHewett et alJ 119951 ) out of 10 4 quasars, and 
among these 16, the confirmed phy sical associations are 
less than 10 ([Kochanc k et al.lll999f) . Even the most re- 
cent optical samples imply the same fraction of binaries, 
0.1-0.4%. Hennawi et al. (2006) analyze a sample of 
~ 59.000 quasars in the Sloan survey, identifying only 26 
close (< 50 kpc) systems. On the other hand, about 30% 
of X-ray detected sub-mm galaxies at z = 2 are in pairs 
(Alexander et al 2005). 

There are only a few well studied cases of ongoing 
mergers at sub-galactic scales in which each companion 
is an AG N; e.g., NGC 6240 (Komossa et al. 2003) and 
Arp 229 (Bal lo et al.ll2004f ). Both are advanced state 
merging systems, belonging to the class of ultraluminous 
infrared galaxies. Noticeably, none of the nuclei show 
signs of activity in the optical band, due to obscuration. 
LBQS 0103-2753 is an example of an optically selected 
double quasar with separation below galactic scale, 3.5 
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kpc (|Junkkarinen et aLll2001[ ). 

Kochanek et al. (1999) analyzed the then available 
optical sample of double quasars in order to assess the 
correspondence between mergers and quasar activity. 
Th ey concluded that the merger scenario (as proposed 
bv lBahcall eTaT][T99l was in good agreement with the 
double quasar population. They supported the merger 
scenario by noticing that the maximum separations be- 
tween pairs are characteristic of scales on which tidal 
perturbations become important (~ 70 kpc), and that 
the absence of small separation pairs is characteristic of 
dynamically evolving binaries. Mortlock et al. (1999) 
analyzes the same sample, providing a dynamical model 
that suggests that double quasars are associated with the 
very first stages of galaxy mergers. They compare the 
dynamical time-scale of the mergers, which results close 
to a Hubble time, to the expected quasar lifetime. Com- 
parison of the time-scales implies that binary quasars are 
short-lived. They further suggest that the most likely ex- 
planation for the death of the binaries is that the in-flow 
of gas ceases, probably as the merger becomes more sta- 
ble. 

We build-up on the inspirational paper by Mortlock 
et al. (1999), and analyze the larger sample of dou- 
ble quasars available today. We focus on the currently 
favored merger scenario, comparing the double quasar 
population to models related to merger simulations, and 
deriving black hole and galaxy properties using the most 
recent results on quasar-host correlations. The redshift- 
dependencies of the observationally-dctcrmincd relations 
linking MBHs to their hosts is an important factor in de- 
termining the merging timescale. We discuss the implica- 
tions of the claimed observed redshift-dependent scaling 
relations in terms of quasar host and pair fraction evolu- 
tion. 

2. SAMPLE SELECTION 

For our sample, we select a to tal of 85 quasa r pairs 
reported in the follo wing pa pers: iHennawi et al.1 (120061) 
H06), IMvers et al.l (|2007a[ ) (M07a). and iMvers et all 
2007bf ) (M07b). The quasars in our sample are selected 
as physical pairs rather than lenses. The authors of our 
sample papers list several ways to establish whether a 
quasar pair is a lens or a binary. Characteristics of a lens 
include similar optical and radio flux ratios between the 
two images, measured velocity differences < 1000km s _1 , 
an observation of the lens candidate, an observation of 
more than two images, a measured time delay between 
images, and an identical set of spectra. 

H06, from which we take 69 out of our 85 quasar 
pairs, report quasar pairs obse rved using three samp les 
from the Third Data Release ([Schneider et al.l 120051 ) of 
the Sloan Digital Sky Survey, namely the SDSS Spec- 
troscopic Quasar Sample, the SDSS Faint Photometric 
Quasar Sample, and the SDSS + 2QZ Quasar Sample. 
H06 select quasar pairs using three different methods. 
For pairs at separation A8 < 3", the authors examine 
the SDSS in search of gravitationally lensed quasars and 
select pairs that do not meet the lens criteria discussed in 
IKochanek e~a l. (1999). For pairs with A9 > 3", the Au- 
thors further compare the emission of the two quasars 
at radio wavelengths. A pair must be a binary if it is 
composed by one radio bright image and one radio quiet 
image. Three of the reported pairs are confirmed as bi- 



naries using this method. For pairs at larger separations, 
they use color selection via the digital photometry of the 
SDSS. Follow-up observations were made to confirm the 
binary hypothesis. 

Finally, M07a report quasar pairs found in the Fourth 
Data Release of the SDSS. T he authors use the co lor 
selection technique outlined in iRichards et alj (|2004t ) to 
distinguish quasars among ~ 300, 000 SDSS objects. 
M07a reports 111 quasar pairs with separations A8 < 
0.1'. Of these 111 pairs, we select eight, requiring both 
images be at the same photometric redshift. We use all 
pairs from Table lb in M07b. These quasars are con- 
firmed binaries and are all included in M07a. 

Table 1 lists the data for each pair. The distribution 
of physica l 3-d separations, assuming a concordance cos- 
mology (jSpergel et al.l [2007) , is shown in Figure [1] We 
do not find any trend between quasar separation and 
luminosity. To estimate the significance of the distri- 
bution over chance superpositions, we subtract in every 
bin the number of companions expected within that bin, 
based on the extrapolation of the quasar correlation func- 
tion. For simplicity, we adopt a single p ower-law func- 
tion £(r) = (r/ro) 7 ([Porciani et al. 2004J, where ro=6.6 
Mpc, 7=- 1. 53), although our sam ple covers a large red- 
shift range (|Hennawi et all 120061 see also). Most of the 
pairs are found either at galactic separations (< 60 kpc), 
in agreement with the results, and the arguments, of 
Kochanek (1999) and Mortlock (1999). The scales over 
which binarity is common correspond to those where 
tidal perturbations become important during the early 
phases of galaxy mergers. A number of pairs above the 
expected number of chance superpositions appears also 
at larger scales (~ 150 kpc). We will discuss these pairs 
further in sections 6 and 7. 

3. SAMPLE ANALYSIS 
3.1. Bolometric Luminosities 

Our first step consists in creating a uniform photo- 
metric sample, by deriving £>-band luminosities for all 
quasars in our sample. 

We convert the g- and u-bands to the £?-band using 
a conversion appropriate to quasar spectra l energy dis- 
tribu tions: B = g + 0.17(u - g) + 0.11 ([Jester et al.l 
2005). From the apparent magnitude B, we derive 
the i ntrinsic luminosity Lb in a concordance cosmol- 
ogy (iSpereel et all 120071 ). The conversion proposed by 
Uester et al.l ( 20051 1 was found using synthetic photom- 
etry using quasars from the SDSS as well as from the 
Bright Quasar Survey. 

Bolometric luminosities L^oi are found using the rela- 
tion: 

logLs = 0.80 - 0.067/: + 0.017£ 2 - 0.0023£ 3 , (1) 

where C = log(L bol /L Q ) - 12 ([Marconi et al.l 120041 ). 
Bolometric luminosities range from 6.69 x 10 n L Q to 
5.94 x 10 14 L Q . The average bolometric correction found 
from the above equation is Li, i/Lb = 5.6. 

3.2. Black Hole Masses 

For a sample of 41 quasars belonging to the Sloan 
Da ta Release 3 we can u se MBH masses derived 
by iVestergaard et all (|2008[ ) by adopting the scaling 
relations based on the broad emission line widths 



3 



2 
I 



25 
20 
15 
10 

5 



20 
10 h 



i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — r 




I i i i i 1 i i i i I 




100 200 300 400 
separation (kpc) 



500 



Fig. 1. — Distribution of physical separations of double quasars. 
The overall distribution peaks at ~ 30 kpc, with a tail extending to 
more than 100 kpc. The bottom panel shows the distribution when 
the extrapolation of the power-law quasar correlation function is 
subtracted (Porciani et al. 2004). Both the peak ~ 30 kpc, and 
the large separation pairs (100 — 200 pc) are above the expected 
number of chance superpositions. 

and nuclear continuum luminosities (e.g. IVestergaard! 



2002t IWarner et al.N2003t IVestergaard fe Petersoni l2006; 
Peterson et al H2001 T For the remaining quasars, we de- 



rive black hole masses by applying a constant Edding- 
ton rati o, fEdd- We comput e fEdd on the sample com- 
piled bv IWoo fe Urrvl (|2002[), selecting sources with lu- 
minosities in the same range as our quasars. We find 
a media n Eddington rat i o /em of 0.29 with a 0.3 dex 
scatter. iKollmeier et~aTl ((2006) find f Edd = 0.25 with 
a 0.3 dex scatter on a sample of X-ray and mid- infrared 
(24/zm) AGNs from the AGN and Galaxy Evolution Sur- 
vey where black hole masses are measured using H/3, 
Mgll, and CIV line emission along with continuum lumi- 
nosities. This result is consistent with our Eddington ra- 
tio of 0.29. Using L bo i/L Q = 3.27 x IQ A f Edd {M BH 
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black hole masses range from 
6.26 ±^3° xl0 10 M o . 

As a further check on our bl ack hole ma s ses, w e apply 
also the relation suggested bv iPeng et alj (|2006l ). based 
on the blue luminosity of the quasar: 



\ogM BH ,% = 7.88 + 0.79 x log 
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where Mbh,8 is the black hole mass in units of 10 8 M C . 



We convert L h „i to ALa.rioo using n hnl ~ »a^ A;5100 
((Peterson et alJl2004 iKaspi et alJl2000f >. The black hole 
mass range is 1.7 ± 1.0 x 1O 8 M to 3.7 ± 3.4 x 10 10 M Q . 
These black hole masses are consistent within error with 
those found from the Eddington ratio. Finally, as a fi- 
na l check, for the 41 black h oles with masses determined 
bv IVestergaard et~a l. ( 2008), the masses we derive step- 
ping through bolometric correction and Eddington ratio 
are within a factor of two from the more reliable masses 
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estimated by IVestergaard et al] (120081). When e ver a n 
MBH mass is available from IVestergaard et~aTl pool) , 
we use the spectroscopically determined MBH value for 
our analysis. 

3.3. Host Galaxy Properties 

Since our goal is to investigate the dynamical prop- 
erties of quasar pairs, we need to embed quasars and 
their black hole engines in appropriate hosts, providing 
the potential for the dynamical decay. Since measuring 
quasar host properties rather than with Hubble Space 
Telescope is daunting, we adopt here an indirect method, 
in which we assume that bla ck hole masses scale with 
the properties of thei r hosts ([Ferrarese fe Merrittl 120001: 
Tremaine etaK 120021: iTreu et al] 12004 iMcLure et all 
20061: IPeng et alj 120061 : IHopkins et alj 120071 ). Albeit in- 
direct, our method allows us to consider how different 
correlations, and their possible redshift evolution, affect 
the dynamical evolution of quasar pairs. 

Using black hole masses, Mbh, we then proceed to 
calculate several properties of the galaxies that host the 
quasars in our sample, namely stellar velocity disper- 
sion (7, stellar bulge mass Mbuige, dark matter halo mass 
Mhaio-, and dark matter halo radius Rhaio- Our starting 
point for the analysis will be in o ne case the Mbh~o~ rela- 
tion fromfTrcmain e~e"t al] (|2002l . case A), a nd in a second 
instan ce the MBH-Mb u i ge relation from lHaring fc Rixl 
(2004, case A), which docs not include redshift evolu- 
tion. We also test the redshift ev olution of the cor rela- 
tion between Mbh~o~ sugg ested bvlTreu et alj (120041) . the 
MBH-Mt, u ig e relation from lMcLure et alj ( 20061 ): cases C 
and D. We will comment on this test in section 4. 

3.3.1. Case A: the velocity dispersion 

Our first endeavor at estimating the properties 
of the host stems from the correlation between 
black hole mass and velocity dispersion of the host 
bulge (iFerrarese fc Merrittl l2000t ITremaine et all 120021 : 
IWoo et al.1 120081 ). In particular, we ado pt the correla- 
tion proposed bv ITremaine et a l. (2002), for the more 
simplified definitio n of veloci t y disp ersion, and the rela- 
tion proposed bv lWoo et al] (|2008l ) to test evolutionary 
effects (see Section 4). 

From black hole masses we compute stellar velocity 
dispersions via: 



logilWg = 0.13 + 4.02 log a 200 . 



(3) 



From this relation, we find stellar velocity dispersions to 
range from 160km s" 1 to 540km s _1 . 

Since one of our goals is to compare the sample of 
quasar pairs to the merger scenario, we relate the bulge 
mass to a velocity dispersion by assumi ng the "funda- 
menta l plane of black holes" proposed by IHopkins et al.1 
(2007), where the authors use galactic merger simula- 
tions to relate MBH masses, host bulge mass and veloc ity 
dispersion. The plane defined in IHopkins et al] (|2007l ) is 
claimed to be redshift-independent, if quasar progenitors 
at each redshift have a wide range of gas fractions. Once 
we derive velocity disper sions from the MBH m asses us- 
ing the relationship from ITremaine et al.1 (|2002l ) we use: 

log M BH ,s = 0.11 + 0.54 log A/ 6u i ffe ,n + 1.82 log er 20 o (4) 

to find the bulge masses of the host galaxies in our sam- 
ple, where 0200 is the velocity dispersion in units of 
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200 kms 1 and Mt, u i g e.ii is the bulge mass in units of 
10 n M o . 

From this relation we find Mb u i ge to range from 3.2 x 
10 10 M Q to 1.2 x 10 13 M o . Assuming that the stel- 
lar velocity dispers ion is a proxy for the virial veloc- 
ity (jFerraresd 120021 ). we can find the properties of the 
host dark matter halos. From the virial theorem and 
the spherical collapse model, we derive dark matter halo 
masses that range from 4.5 x 10 12 M Q to 1.9 x 1O 13 M , 
and the radii of the dark matter halo are in the range 
190 kpc to 660 kpc. The masses of the host halos are in 
agreement with, although slightly higher than, the typi- 
cal masses of quasar hosts derived from clustering anal- 
ysis (jMvers et al.ll2007at ICoil et al.ll2007l see Fig. 2). 
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Fig. 2. — Physical separations of double quasars as a function 
of redshift (solid circles). The sizes of the halos that we derive 
for our sample are shown as well. Open circles: case A, velocity 
despersion. Asteriscs: case B, bulge masses. For comparison we 
show by a solid curve the size of a halo with mass 5 X 10 Mq, the 
typica l quasar host hal o found in clustering analysis IJVlvers et alj 
2007a; Coi l et al.|[2007Tl . 



3.3.2. Case B: the bulge mass 

We can derive the mass of the host bulge by apply- 
ing two observational, empirical relations, one that was 
found for loca l quie scent MBHs (|Marconi fc HuntJ 120031: 
Marine fe Rixl l200l and one that was found for high- 
redshift radio-loud quasars iMcLure et al.l (120061). Our 
case B adopts the relationship proposed bv lHaring fc Rixl 
(HoU): 



\ogM BH ,8 = 0.2 + 1.12 log M bu . 



(5) 



Bulge masses are in the range 2.9 x 10 Mq to 6.2 x 
10 12 M Q . We will disc uss the effect o f adop ting the re- 
lationship proposed by IMcLure et al.1 (|2006f ) in the next 
section. 

We now reverse the method used for case A and we ap- 
ply Equation 2 to compute velocity dispersions starting 



from the stellar bulge masses. Velocity dispersions range 
from 170km s -1 to 840km s _1 . We can compare the ve- 
locity dispersion that we d erive "directly" from th e black 
hole properties, e.g., using iTremaine etafl (|2002f ). to the 
velocity dispersion that we deri ve "indirectly" by app ly- 
ing in order the relationship by Haring fc Rixl (|2004l ) to 
find the bulge mass and then the Black Hole Fundamen- 
tal Plane to find the velocity dispersion, from black hole 
and bulge mass. The relationship reads log Mbh, 8 — 
0.02 + 3.50 log 0200, leading to a steeper dependence of 
a on BH mass. The "indirect" method therefore re- 
sults in a larger velocity dispersion at fixed BH mass 
for Mbh > 2 x 10 7 . From the stellar velocity dispersion 
we once again find the dark matter halo mass, ranging 
from 4.8 x 10 12 M o to 3.0 x 1O 13 M . The radii of the 
dark matter halo range from 180 kpc to 840 kpc. 

4. REDSHIFT EVOLUTION OF MBH-HOST 
CORRELATIONS 

In the previous models we have assumed that the cor- 
relations between black holes and their host are redshift- 
independent. This is indeed clai med, theoretically, fo r 
the black hole fundamental plane (jHopkins et al.1 1 2005). 
However, the two correlations between Mbh and M{, u i g e, 
and Mbh and a might evolve with the cosmic time 
(IMcLure et al.1 120061 ; iPeng et al.|[200l iTreu et al.112001 
IWoo et al.H2006D . 

4.1. Case C: redshift evolution of Mbh — Mb u i ge 

IMcLure et all (|2006l) use the 3C RR sample of radio- 
loud quasars to investigate the evolution of the correla- 
tion between black hole and bulge masses in the redshift 
range < z < 2. Stellar bulge masses are computed by 
the authors for a sample of radio-loud quasars from the 
3C RR sample from if-band magnitude measurements, 
while black hole masses are found from H/3, Mgll, or CIV 
line widths. Mb u i ge -z and Mbh-z relations are found in- 
dependently, and the two are combined to form the rela- 
tion we use. Below z = 1 the 3C RR sample is consistent 
with a black hole weighing 0.002 times the mass of the 
bulge. The ratio between black hole and bulge mass is 
found to evolve with redshift as Mbh /Mb u i ge ~ (1 + z) 2 . 
We define a case C, where from black hole masses we 
compute stellar bulge masses via 

\og{M BH /M bu i ge ) = 2.07 log(l +z)- 3.09 (6) 

(jMcLure et al.ll2006f ) , and subsequently calculate velocity 
dispersions using Equation 4. 

4.2. Case D: redshift evolution of both Mbh — Mbuige 
and Mbh — & 

ITreu et al.1 (|2004l ) and lWoo et all (|2006D estimate black 
hole masses for a sample of AGN at z = 0.36 from the Hp 
line width and the optical luminosity at 5100 A, based on 
the empirically ca librated photoionization method (e.g., 
iKaspi et al.l l2000h . They also directly measure velocity 
disp ersions from stella r absorption lines in the host spec- 
tra. IWoo et all (|2006D find that, at a fixed velocity dis- 
persion, the mass of the MBH increases with redshift. In 
other words, MBHs are hosted by galaxies with smaller 
velocity dispersion at earlier cosmic times. At redshift z, 

\ogM BH {z) - logAfsff (0) = 1.66* + 0.04. (7) 



Fig. 3. — Comparison between galaxy properties derived directly 
from black hole masses, using the relationships found for z > 
black holes (Woo et al. 2008; McLure et al. 2006). Left panel: 
velocity dispersions. Right bottom panel: bulge masses. Right top 
panel: bulge mass vs halo mass, if we assume that the velocity dis- 
persion can be used as a tracer of the total potential (e.g., Ferrarese 
2002). 

iTreu et all ()2007f ) further find that, for the same galax- 
ies in the sample, MBHs are hosted by galaxies with 
small er bulges, with a r edshift dependence similar to 
what lMcLure et all (|2006D fi nd (at least for z < 1). The 
redshift scaling proposed bv lWoo et alj (|2006f ) has been 
derived on a sample at a very specific redshift, much 
lower than the redshift of our quasars (see Table 1 and 
Figure 2). The extrapolation of their relationship at 
much higher redshifts leads to some "peculiar" cases, 
such as our highest redshift quasar, apparently powered 
by a 1.4 x 10 10 M Q MBH, strikes us as being hosted by a 
bulge with velocity dispersion less than 50kms _1 . Fur- 
thermore, in this case if we were to assume that the veloc- 
ity dispersion can be used as a tracer of the dark matter 
potential the mass of the halo would be roughly equiva- 
lent to the mass of the bulge. 

A more recent analysis by I Woo et a l. (2008) modifies 
the redshift scaling as: 

log M BH (z) - logAfsff (0) = 3.1 log(l + z) + 0.05. (8) 

This scaling relation improves somewhat on the "pe- 
culiar" cases mentioned above (now the host of the 
1.4 x 10 10 M MBH has a velocity dispersion of ~210 km 
s _1 ), and the galaxy properties are now illustrated in Fig- 
ure [3] Still, the redshift dependence leads to very differ- 
ent properties for the hosts, which reflect on the dynam- 
ical evolution, when we compute the merging timescales. 
Figure [3] illustrates these differences. If the relationships 
between black hole masses and their hosts are redshift- 
dependent, typically MBHs are hosted in smaller galaxies 
(both smaller velocity dispersions, and smaller bulges). 
In most cases, therefore, the separation between the 
quasars in a pair are larger than the size of their host 
halos, leading to the conclusion that the merger of the 



Fig. 4. — Dynamical friction timescale for quasar pairs, in case A 
(bottom) and case B (top). Solid curves are linear fits. The dotted 
horizontal line marks the Salpeter time, tsalp ~ 4.5 X 10 7 yr(e/0.1). 

systems is still in its first phases, and has a long way 
until coalescence. 

Our cas e D couples the rel ationship by Woo et al. 
(2008) and lMcLure et alj (|2006f ) to find joint velocity dis- 
persions and bulge masses. 

5. DYNAMICAL FRICTION TIMESCALES 

From the host galaxy properties, we can compute dy- 
namical friction timescales t^F, i.e. the time needed for 
the host galaxies to merge, using the Chandrase khar for- 
mula 2 (M99, see also lBinnev fc Tremaine]|1987l ) 

t DF ~4x lO 9 (r ma;c /2Okpc) 2 (a 2 oo)/(lO 9 M /M)yr. (9) 

This equation is derived by integrating the equation 
for the time spent at a separation r, \dt/dr\ ~ 
3.3(a/G?Ml n(A))r, where In (A) ~ ln(2r»er 2 /GM) and 
r» ^ lOkpc (ILacev fc Cold lT993). Here a is the velocity 
dispersion of the largest host galaxy, and M is the mass 
of the infalling system. To find tjjp for the galactic merg- 
ers of our sample, we use the separation of the quasars 
as r max and the largest velocity dispersion in the pair 
as a. M has to be related to properties of the smaller 
system (the "satellite"). For quasars with separations 
larger than the size of their host halos, we assume that 
the merger of the systems is still in its first phases, and 
we expect the total time to coalescence and formation of 
the new galaxy to be determined by the global potential 
wells, hence we choose M = Mhaio in equation [5J When 

2 Detailed calculations of the merging timescales will depend on 
the galaxy potent i als an d on orbital parameters of the encounters, 
e.g. ITaffoni et all pOOST) : IBovlan-Kolchin et alj gOOS). Given the 
uncertainties on which is the most appropriate choice for the poten- 
tial wells of galaxies, and the fact that there is at most a factor of a 
few diskrepancy between the dynamical friction timescale and the 
merging timescale in the case of major mergers, we prefer a sim- 
ple, clear, physically motivated choice such as the Chandrasekhar 
formula. 
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Fig. 5. — Dynamical friction timescale for quasar pairs, in case 
C (bottom) and case D (top). 

instead the two quasars appear to be already "inside" 
the main host halo, we use M = M bu i ge , assuming that 
the blac k holes ar e still surrounded by their stellar en- 
velopes (|Yu|[200l . Figure[2] compares the separations 
between pairs to the typical halo sizes for cases A, B. 
Only about 20% of mergers appear to be still "outside" 
the respective halos of galaxies. 

Our treatment is very simplified, as it does not 
include the possibility that we are detecting the 
quasars at pericenter or a pocenter on eccentric orbits 
(jKhochfar fc BurkertJ 1200a ) . however, it is very unlikely 
to detect a merger at pericentric passage, where the sys- 
tem spends the least amount of time, and our results are 
not strongly biased if the systems are near apocenter. 
The timescales we find for galactic mergers are in rough 
agreement with simulations when we apply the appropri- 
ate radii and masses from our sample. For both cases A 
and B we find that the timescale s range from 10 6 yr t o 
10 10 yr, compared to ~ 10 9 yr from lHopkins et aT] (2005). 
The agreement is best at low redshifts, with merger 
timescales decreasing with increasing redshift both in 
case A and B. We show the merger timescales of the 
quasar hosts as a function of redshift in Figure [4] In 
our cases C and D the merger timescales vary extremely 
slowly with cosmic time (Figure [5]). The typical merging 
timescales remain of the order of 10 8 — 10 10 



years. 



6. RESULTS 



We compare the merger timescales to the lifetime ex- 
pectancy of quasars in order to understand the link be- 
tween mergers and quasar fueling: N dmMe ^ qso /N qso ~ 
tqso/tDF (cfr. Martini 2004). A numerical investigation 
of the qu asar lifetime in galax y mergers has been pre- 
sented in Ho pkins et al.l (|2005D . They suggest that the 
time spent at a given luminosity is "luminosity depen- 
dent". The quasars in our sample have Lb > lO 11 !/©, 
which leads to t qso < 10 7 yr. However, luminosities larger 
than Lb > 10 11 L Q are predicted to occur only after the 



systems have coalesced, in the merger remnant, where 
the quasar separation is well below the SDSS resolution 
limit. Since this set of simulations is not representative 
of the double quasar population that we are investigat- 
ing, for sake of simplicity we adopt as "standard" value 
associated to quasar lifetimes the Salpeter time, is a i P = 
eccT / \ATrGm p ) ~ 4.5 x 10 7 yr(e/0.1), where e is the radia- 
tive efficiency. If we compare t qso ~ 4.5 x 10 7 yr to the av- 
erage merging time scale, we find that N douU e,qso/N qso ~ 
4 — 5% for case A ((tnj?) = 9.2 x 10 8 yr), and similarly 
for case B ((t DF ) = 1.3 x 10 9 yr). 

Figure 2] also suggests that, depending on what the 
intimate link between black holes and their host really 
is, we should see a varying fraction of double quasars 
with redshift. We can compute the average merg- 
ing timescale in 3 redshift bins, and repeat the same 
analysis. In case A, for quasars at z < 1 we find 
(tDF) = (4-1 ± 6.2) x 10 9 yr, yielding a binary fraction 
NdouUe,qso/N qso < 1%. For quasars at 1 < z < 2 the 
fraction increases, as (toF) = (3.5 ± 5.6) x 10 8 yr, yield- 
ing a binary fraction N double, qso / 'N qso ~ 3%. In our last 
redshift bin, 2 < z < 3, the average merging timescale is 
similar to the previous bin, but the dispersion decreases, 
as more quasars have similar, short, merging timescales: 
(t DF ) = (3.2 ± 3.9) x 10 8 yr, so N double , qso / N qso ~ 3%. 
Results are similar for case B. In our cases C and D 
the merger timescales are roughly redshift-independent, 
leading to the expectation of an approximatively con- 
stant binary fraction 

We postpone a detailed analysis of the redshift depen- 
dence of quasar pairs to a forthcoming paper, for now 
we show in Figure [5] some tentative evidence that the 
pair fraction might indeed be redshift dependent. Here 
we can calculate the binary fraction in the same redshift 
bins where we have compared t qso to the merging time 
scale. We find N doub i e ,qso/N qso < 0.1% at z < 1, and 
also at 1 < z < 2, while N double ,q S o/N qso ~ 0.2% at 
2 < z < 3. 

Assuming that only one physical process triggers 
quasars, we can compare the overall redshift distribu- 
tion of quasars in a homogenous statistical sample to 
the distribution of pairs only. We select only double 
quasars from Hennawi et al. (2006) and compare their 
redshift distribution to that of all quasars in the Fifth 
Data R elease of the Sloan D igital Sky Survey quasar 
catalog ([Schneider et al.l I2005T ). Figure [5] shows a hint 
of increasing pair occurrence at z > 1.5. 

Finally, we compare the typical separations at which 
binary quasars are observed (Figure [1]) to the merger 
phase and separations that should correspond to the 
quasar phase in merger simulations. We find that the 
distribution of physical separations peaks at ~ 30 kpc, 
with a tail at larger separations (~ 100 — 200 kpc). The 
peak of the distribution is roughly consistent with the 
first, earliest (low-luminosity) episode of quasar activity 
found in equal mass mergers simulations (see Figure 2 in 
Hopk ins et al .1120061 ) , although the numerically predicted 
luminosities in this phase of the merger are smaller than 
the typical double quasar luminosity in our sample. The 
tail of the quasar pairs distribution at large separations 
is instead inconsistent with any quasar activity predicted 
by published simulation results. 

Nevertheless, the mass ratio of the black holes in the 
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Fig. 6. — Redshift distribution of double quasars in our sam- 
ple. Bottom panel: all quasars in the sample. Top panel: double 
quasars from Hennawi et al. (2006) only (solid histogram). For 
comparison we show the redshift distribution of all quasars in the 
Fifth Data Release of the Sloan Digital Sky Survey quasar cata- 
log (Schneider et al. 2007, dashed histogram), renormalized by a 
factor of 10 -3 . 

pairs (that we can use as a proxy for the typical mass 
ratios of the merging galaxies) does not match the 1:1 
ratio that was typically selected in galaxy merger sim- 
ulations. The average mass ratio is (2.4 ± 1.3), which 
corresponds to a regrettably unexplored range. When 
we consider a smaller secondary galaxy, we expect a 
smaller deceleration due to dynamical friction, and con- 
seq uently a larg e r gala xy separation at the first apocen- 
ter. iColpi et alJ (l2007[) indeed sh ow that in a 1:4 galaxy 
merger ([Kazantzidis et all 120051 ) the first apocenter oc- 
curs at approximately 180 kpc. 

The presence of the large separations tail suggests that 
a significant fraction of double quasars is seen at the first 
apocentric passage, where the two galaxies spend a much 
longer time with respect to the pericenter. This scenario 
would follow naturally if gas is dynamically perturbed 
during the first pericentric passage, forming soon after 
a gaseous circumn uclear disk (with a size of a few hun- 
dred parsecs, e.g. Kazantzidis et al.l 120051 ), but a delay 
occurs before the gas in the disk looses enough angular 
momentum to be driven at subparsec scales where it can 
finally be accreted by the MBH. This final stage likely 
occurs during the first apocenter, where the quasars are 
observed, causing the tail at ~ 100 — 200 kpc. 

7. CONCLUSIONS 

If most galaxies host a MBH (e.g., Richstonc 2004) and 
galaxy mergers trigger quasar activity (Haehnelt & Rees 
1993; Wyithe & Loeb 2003; Kauffmann & Haehnelt 2000; 
Cavaliere & Vittorini 2000; Di Matteo et al. 2005) then 
the observational sample of binary quasars can be used 
as a probe of quasar lifetimes. Let's start with Occam's 
razor and assume that every galaxy hosts a MBH and 
quasar activity is triggered by galaxy mergers. If the 



lifetime of quasars equals the merger timescale, t qso = 
toF, the probability of observing a double quasar is in 
principle 100%, if we do not consider additional factors, 
such as obscuration. If t qso « tpf , then the probability 
much smaller than 100%. Additionally, if there is a delay 
in the triggering of the two quasars, then one might have 
ceased its activity before the other started, for instance. 
Face value, our results require that the lifetime of quasars 
is extremely short, about an order of magnitude shorter 
than the Salpeter time, or the lifetime expected from 
numerical simulations. 

In a merger driven scenario for quasar activity, with 
MBHs in most galaxies, where naively all sources are 
expected to be in binaries, the 0.1% binary fraction can 
be explained in different ways: 

• Quasar pairs are obscured much more than ex- 
pected from simulations, at Compton thick levels 
(Njj > 10 cm -2 ), otherwise pairs would be easily 
detected in X-rays. 

• Activity is triggered on only one of the 2 black 
holes, or on the MBH(s) in the merger remnant. 
Even if there is a MBH binary in the remnant (Hop- 
kins et al. 2006), we should expect the binary to be 
hard, or at least the MBHs to be bound by the time 
the galaxy merger remnant has relaxed enough to 
be considered a " normal galaxy" . A pc-scale binary 
would not be resolved in the SLOAN or LBQS sur- 
veys. However, this hypothesis does not explain 
pairs occuring at ten-kpc or hundred-kpc scales. 

• There is a time delay between the onset of activ- 
ity in the 2 merging systems (e.g. the gas infall 
timescale depends on the dynamical time of each 
host, rather tha n the dynam i cal tim e of the merg- 
ing system, see IColpi et all (|2007l )L If the two 
quasars light up at different, random times, then 
N d oubie,qso/N qso ~ (t qso /t DF ) 2 ~ 0.1%, consistent 
with the observational result. If the time delay is 
longer that the quasar timescale, then binarity is 
avoided completely. 

The average mass ratio of the quasar pairs we analyzed 
is ~ 2, while, since the mass function decreases steeply at 
large galaxy masses it is more probable that a very mas- 
sive galaxy 3 interacts with a much smaller secondary. If 
the mass ratio is very large, quasar activity is not sup- 
posed to be triggered at the level of observability our 
sample requ ires. If the mass r atio is > 4, quasar activity 
is expected (jColpi et al.1120071 ). but the delay between the 
onset in the primary and in the secondary is too long for 
the pair to be detected as a "double quasar" . Hence we 
prefer the third hypothesis above. 

Clearly the occupation fraction (that is the fraction of 
galaxies hosting a MBH) is indeed an important factor. 
If not all galaxies host a MBH the fraction of binaries 
decreases. As long as we do not have a demography of 
inactive MBHs at z > 0, this issue cannot be addressed 
properly. This fair comparison has to await for the Laser 
Interferometric Space Antenna, probing the BH popula- 
tion via gravitational waves detections. However, we can 

3 As all the MBHs in our sample have masses above 10 8 Mg , 
they are hosted in large galaxies sitting at the steep end of the 
mass function. 



compare our hypothesis to theoretical expectations. In 
quasar evolution models that evolve the population ab 
initio a minimum number density of MBHs is required 
to match the luminosity function of quasars. That is, 
one can not have less MBHs than quasars, at a given 
time. With this constraint in mind models such as those 
described in Volonteri et al. 2008 (see also Marulli et 
al. 2007) require an occupation fraction higher than 
90% for MBHs hosted in halos with mass above 10 12 M Q , 
the typical quasar host halo found in clustering analysis 
(jMvers et al.ll2007aHCoil et al.|[2007l ). 

We compared the typical separations at which binary 
quasars are observed (Figure [T]) to the merger phase and 
separations that should correspond to the quasar phase 
in merger simulations. We find that the distribution of 
physical separations peaks at ~ 30 kpc, with a tail at 
larger separations (~ 100 — 200 kpc). The peak of the 
distribution is roughly consistent with the first episode 
of quasar activity fou nd in equal mass me rgers simula- 
tions (see, Figure 2 in [Hopkins ct al. 2006|). The tail of 
the quasar pairs distribution at large separations does 
not correspond to any quasar phase predicted by pub- 
lished simulation results. However, the mass ratio of the 
black holes in the pairs does not match the 1:1 ratio that 
was typically selected in galaxy merger simulations. The 
average mass ratio is (2.4 ± 1.3), a regrettably still unex- 
plored range. 

The pairs at large separation could be unequal mass 



mergers seen at the first apocentric passage. This sce- 
nario would follow naturally if gas is dynamically per- 
turbed during the first pericentric passage, but the gas 
in the disk requires longer timescales to loose enough an- 
gular momentum to be driven at subparsec scales and 
fuel the MBH. This final stage likely occurs when the 
two galaxies are at the first apocenter. 

Finally, we explored how the redhift dependence (or in- 
dependence) of the correlations between MBHs and their 
hosts can be tested using large samples of quasar pairs. 
If MBHs follow a redshift independent correlation with 
either the velocity dispersion of the host (e.g., Tremaine 
et al. 2002, Ferrarese & Merritt 2000, case A), or with 
the mass of the bulge (e.g. Marconi & Hunt 2004, Har- 
ing & Rix 2004, case B), we expect merger timescales to 
decrease at high redshift, leading to a higher fraction of 
double quasars. If, on the other hand, the correlations 
between MBHs and their hosts evolve with time (cases 
C and D, e.g., McLure et al. 2006, Woo et al. 2006) 
the merger timescales vary extremely slowly with cos- 
mic time, and an increase in the binary fraction is not 
expected. 



We wish to thank Marianne Vestergaard for kindly 
providing black hole masses for SDSS Data Release 3 
quasars. 
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TABLE 1 

Properties of the double quasar sample. We list the name and 
redshift of the pair, the physical separation between the two 
quasars, along with the luminosity and bh mass for both quasars in 
the pair. the last column lists the reference for the photometric 
or spectroscopic information. (1) hennawi et al. (2006), (2) myers et 
al. (2007a), (3) Myers et al. (2007b), (A) Vestergaard et al. (2008) 
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T f O 1 A AA 

1.68 ±0.09 


bDbbJ0054-0946 


2 


13 


117.11 


2.63+0. 151j+1 


A AA 1 A T T 
2.02 ±0.11 


O T^i OD TA1 1 ^7 i n ro 

bDbbJ 01 17+3153 


2 


13 


93.85 


2.38 ±0.13 


A A O 1 A A P" 

4.43 ±0.25 


O T"\ O O T a a a i i aaoa 

bDbbJ 0201+0032 


2 


30 


155.90 


7.63 ±0.42 


o o r 1 A T A 

3.35 ±0.19 


t < t \ ( ■ ( < t a a j o i a a a / \ 

bDbbJ 0248+0009 


1 


64 


58.45 


a r\r7 i r\ no 

4.07 ±0.23 


T AT 1 A A/' 1 

1.07 ±0.06 


o T"\ o o Tnooo / \ •-• 1 1 ,- 1 

bDbb J 0332-0722 


2 


10 


150.62 


2.95 ±0.16 


A A A 1 A T A 

2.24 ±0.12 


O T^i OD TAO J/ 1 1 AT A a 

bDbb J 0846+2749 


2 


12 


156.25 


5.84 ±0.32 


A A T 1 A A A 

4.41 ±0.24 


O T"\ O O T/ \(\.>(1 1 P" A P" O 

bDbbJ0939+5953 


2 


53 


262.47 


6.34 ±0.35 


T O F" 1 A T ATT 1 1 T 

1.85±0.10E+1 


bDbbJ0955+6045 





72 


134.40 


1.81±0.101j-1 


T P"A 1 A AOTT 1 T 

1.50±0.08E-1 


o T""\ o o Trifirn i p" a ac\ 

bDbb J 0959+5449 


1 


95 


32.73 


2.93 ±0.16 


T O r' 1 A T A 

1.86 ±0.10 


O T~\ O O T 1 A 1 A l C\ A~\ C 

bDbb J 1010+0416 


1 


51 


145.59 


1.56 ±0.09 


T A T 1 A A O 

1.41 ±0.08 


ot^oo t "i a t a i c\ i \ o r\ 

bDbbJ 1014+0920 


2 


29 


180.65 


9.19 ±0.51 


O O A 1 A T A 

3.39 ±0.19 


onoo t "i aao i o c\ o c\ 

bDbbJ 1028+3929 


1 


89 


63.12 


1 A1 I O 1 O 

3.21 ±0.18 


T OO 1 A AO 

1.38 ±0.08 


o t^i ti ao /i i a ta t 

bDbbJ 1034+0701 


1 


25 


25.85 


1.63 ±0.09 


O OT 1 A T ATT 1 T 

3.37±0.19E-1 


bDbbJ 1056-0059 


2 


13 


59.80 


4.11 ±0.23 


T A A 1 A T T 

1.94 ±0.11 


o t^i o o t t -i a o i a a o t 

bDbbJ 1123+0037 


1 


17 


465.12 


2.22 ±0.12 


{ • C\C 1 A OATT 1 T 

6.96±0.39E-1 


O T~\ O O T1 IOC I tiC A A 

bDbbJ 1225+5o44 


2 


38 


48.92 


o ah n /in 

8.93 ±0.49 


a Ar 1 A T r* 

2.95 ±0.16 


bDbbJ 1254+ol04 


2 


05 


146.91 


O 1 1 A *1 C 

8.33 ±0.46 


PT OA 1 A AA 

5.32 ±0.29 


o t^i T1 orn i "i a a t 

bDbb J 1259+1241 


2 


19 


29.78 


/I (~\f 1 A AO 

4.06 ±0.22 


A AC\ 1 A A A 

4. 42 ±0.24 


O OO T1 1AH i n AA 

bDbb J 1303+5100 


1 


68 


32.18 


1 i ' 1 ~\ 1 A AA 

1.60 ±0.09 


T AA 1 A AT 

1.22 ±0.07 


o T"\ o o t t o t a i z 1 aao 

bDbb J 1310+6208 


2 


06 


391.25 


"1 "1 A 1 A A i~ k l 1 

1.12 ±0.061j+1 


A A T 1 A T O 

2.27 ±0.13 


O O O T T O O T i /^A1 A 

bDbbJ 1337+6012 


1 


73 


26.23 


O TO 1 A JO 

8.73 ±0.48 


AOT ] A T O 

2.31 ±0.13 


OT\0 O T1 O I A i 1 A A*"7 

bDbbJ 1349+1227 


1 


72 


25.39 


1 OA I A 1 ATT 1 1 1 

1.80 ±0.101j+1 


A AC A rip 

4.46 ±0.25 


crxoc n a aa i i aoa 

bDbbJ1400+12o2 


2 


05 


121.87 


O TO A 1 c 

2.73 ±0.15 


A A A 1 A T O 

2.40 ±0.13 


O T^i OO T1 *l AT I A A A T 

bDbbJ 1405+4447 


2 


23 


61.04 


A O ■ 1 A "1 ATT 1 1 1 

2.26 ±0.121j+1 


O /^T 1 A AA 

3.67 ±0.20 


t^i t 1 a r\c\ 1 / \ 1 / \ 

bDbbJ 1409+3919 


2 


08 


56.66 


2.96 ±0.16 


T f T 1 A AA 

1.61 ±0.09 


bDbb J 1530+5304 


1 


53 


34.72 


"1 AA 1 A i \ / ' 

1.02 ±0.06 


A p" r 1 a r oT7 t 

9.55±0.53E-1 


bDbb J 1546+5134 


2 


95 


OA/ 1 / • ( 1 

326.69 


O AT 1 A AA 

3.91 ±0.22 


T AA 1 A ATTT 1 1 "1 

1.22±0.07E+1 


O T^i OO T1 /'AA I 0^70 J 

bDbbJ 1629+3724 





92 


34.49 


A P" f~* 1 A T* O TT 1 T 

9.56 ±0.53b-l 


O A O 1 A A TTT 1 T 

8.48±0.47E-1 


bDbb J 1719+2549 


2 


17 


121.76 


J A^ 1 A A J 

4.26 ±0.24 


A AA 1 A AO 

4.99 ±0.28 


bDbb J 1723+5904 


1 


60 


31.35 


J A A 1 A A r" 

4.49 ±0.25 


1.16 ±0.06 


bDbb J 2128-0617 


2 


07 


69.20 


J AT 1 A AT 

4.91 ±0.27 


A AO 1 A A P" 

4.48 ±0.25 


O T~\ O O TAAT /I 1 "1 O / ■ 

bDbb J 2214+1326 


2 


00 


48.55 


A T/"' 1 A t r 

2.76 ±0.15 


A A A 1 A T A 

2.24 ±0.12 


bDbb J 2220+ 1247 


1 


99 


133.17 


O OA 1 A T O 

3.22 ±0.18 


T PT A 1 A AO 

1.52 ±0.08 


bDbbJ0740+2926 





98 


20.72 


A OC 1 A T O 

2.36 ±0.13 


7.54±0.42.h-l 


o t^i oo n aot i a t p" a 

bDbbJ 1035+0752 


1 


22 


22.45 


A AA 1 A T A 

2.22 ±0.12 


/ • r* -\ I r\ OTTT 1 "1 

6.61±0.37E-1 


bDbb J 1124+5710 


2 


31 


18.04 


T J A 1 A AOT? 1 T 

1.49 ±0.081j+1 


r~ A A | A AA 

5.24 ±0.29 


ot^i 00T1100 i /"*o at 

bDbbJ 1138+6807 





77 


19.26 


T AO 1 A T T 

1.98 ±0.11 


O TP" 1 A AT TT 1 T 

3.75±0.21E-1 


OT^OO T "1 r AC! i OOAO 

bDbbJ 1508+3328 





88 


22.45 


nil 1 A T A 

3.41 ±0.19 


O "1 T 1 A "1 TTT 1 1 

3.11±0.17E-1 


OrvCC TAA1 A i AArri 

bDbb J 0012+0052 


1 


63 


135.55 


OCA 1 A AOT? 1 


A O A I A tr ATT 1 T 

9.o4±0.52.h-l 


OT~\00 TA1 "1 T 1 AAAA 

bDbb J 01 17+0020 





61 


O A A TO 

299.78 


T A A 1 A A O 

1.49 ±0.08 


T T A 1 A "1 ATT 1 T 

1.74±0.10E-1 


OROO TA1 A T I AA11 

bDbbJ0141+0031 


1 


89 


361.05 


A A T 1 A T /I 

2.47 ±0.14 


T /^T 1 A AA 

1.67 ±0.09 


OL/OOJUZ'lO-UllO 


2 


A(\ 
*±v 


^fi A^ 


cc 70 _Lf| ^9 


^ f\A 4-0 1 7 


SDbSJ0258-0003 


1 


32 


246.69 


3.27 ±0.18 


1.40 ±0.08 


SDSSJ0259+0048 





89 


152.21 


8.59±0.48E-1 


1.07 ±0.06 


SDbSJ0350-0031 


2 


00 


380.88 


2.82 ±0.16 


5.66 ±0.31 


SDSSJ0743+2054 


1 


56 


300.71 


2.03 ±0.11 


1.64 ±0.09 


SDSSJ0747+4318 





50 


56.16 


2.17 ±0.12E-1 


1.78±0.10E-1 


SDbSJ0824+2357 





54 


94.66 


5.18 ±0.29E-1 


4.39±0.24E-1 


SDSSJ0856+5111 





54 


138.49 


5.00±0.28E-1 


1.91±0.11E-1 


SDSSJ0909+0002 


1 


87 


126.35 


6.16±0.34E+1 


2.72 ±0.15 


SDSSJ0955+0616 


1 


28 


368.00 


6.99 ±0.39 


8.58±0.48E-1 


SDbSJ1032+0140 


1 


46 


465.77 


4.42 ±0.25 


1.02 ±0.06 


SDbSJ1103+0318 


1 


94 


481.14 


1.43±0.08E+1 


7.58 ±0.42 


SDbSJ1107+0033 


1 


88 


208.81 


7.70 ±0.43 


2.83 ±0.16 


SDbSJ1116+4118 


2 


99 


106.41 


2.63±0.15E+1 


1.12±0.06E+1 


SDbSJ1134+0849 


1 


53 


229.47 


3.90 ±0.22 


2.72 ±0.15 


SDbSJ1146-0124 


1 


98 


238.83 


2.55 ±0.14 


2.15 ±0.12 


SDSSJ1152-0030 





55 


187.86 


4.19±0.23E-1 


1.23±0.07E-1 


SDbSJ1207+0115 





97 


281.41 


1.49 ±0.08 


3.81±0.21E-1 



4.99±§;i 
3.60±1.31 
2.32±' ' 



% 



9-47±i g 46 
2.52±0.55 
2.95±109 
9.87±^# 

1.23±g;ljE+ 

L.y? 



2.24± 



1.53 



6.46±2.53E-1 
1.17±0.77 
2.28±°;| 
2.80±?|? 
1.60±f;| 
1.86±fp 

4.14±0.62E-1 

5.79+1.93E-1 
2.55±0.34 

1.76±|]2iE-l 

l-96±?34 
5.25±^<»E-1 

2.26±S1?E-1 

1A4± fi 

3.27± 2 '24 

8.36±2.13E-1 
1.60±0.94 



1.64±f;i2 



.12 

3.63±1.52E-1 
2.85±Sg^E-l 
3.23±S : }|E-1 
1.79±g;»E-l 
5.92±6;i 
2.08±0.39 
7.14±1.54E-1 

4.32±2 



-1 



i 43^2.r2jr_ 
1.19±0.36 
1.49± 2 ;*f 
5.01±0.88E-1 
2.95±0.57E-1 



1.93± 



7.32±1.16E-1 
5.67±^| 1 E-1 
3.57± 2 -iiE-l 
2 24±» 2 5 

1-39± 2 ;| 
4.28±|-I? 

3.30±*jg 

1.15±0.40 
8.01±4.26E-1 

1.69±?; 2 1 
3.40±2.16E-1 

1.17±0.26 

1.85±fp 

1.43±0!22 
2 ^5+4.48 
z -' 30=l= 1.61 

4.83±1.09E-1 4.27± 



2.23±0.29E-1 
1.18±f|?E-l 
9.93±i 8 7 92 E-l 

1 iq+2.26 
0-81 
2.52±f| 2 

.kit 



Ml 



1.31±?- 



■1 



2.52±f|l 
1.02±J-«E+1 
9.98±I 9 7?°E-2 

1.09±f;0| 
8.37±J 5 7? l:! E-l 

L95 ±L33 
8.21±i 5 J-*E-l 

2.15± J f 



5 4 5 & E . 



1.14# 17 



..78 

4 E- 



1.71+ ' 2: ' 



3.02±? 



2.52±, 
7.31±f| 
1.32±g- 52 

U2 



1.35± 
2.55± 
1.40± 
2.11± 



2 

2 



9.53±i 8 4„ 
5.78±i 1 [ ;S 1 E-l 

fi 77+12.90 

516 7 ±ffil 

2.84±5;| 2 
6.98±J 3 4 ff E-l 

1.4.88 



-Q.89 
¥ E- 



2.56±* 

4.61±|«E-1 
4.07±|:«E-1 

2.97±|^ 
2.37±f jy-E-1 
1.99±? : I§E-1 
**E-1 

1e- 



5.66± 
1.15± 



-1 



9.88±I S 7f s E-l 



1.76±H§ 
8.32±i 5 gf E-l 
6.42±P4 3 E-1 

3.21±f-| 2 
9.70±I 8 f jf E-l 
1.17±|'||E-1 
2.75±? : 11e-1 

1 43± § 54 
9.41±i 5 | 95 E-l 
5.21± 9: ||E-1 
6.14±ii i | b E-l 
4.88±2.23 

i.64±f : « 



6.24± 



1.58±5| 
i 2fi+ 2 " 9 
8.33±| 5 f j| jr E-2 
2.41±j : ||E-l 



1, A 
1 

1, A 
1, A 

1 

1 

1, A 
1 
1 

1, A 
1, A 
1 

1, A 
1, A 
1 

1, A 
1 

1, A 
1 

1, A 
1 
1 

1, A 
1, A 

1 

1 

1 

1 

1 

1 

1, A 
1, A 
1, A 

1 

1 

1 

1, A 
1, A 
1, A 
1, A 

1 

1 

1, A 
1 
1 

1, A 
1, A 
1 

1, A 
1, A 
1 

1, A 
1 

1, A 
1, A 
1, A 

1 

1 

1, A 
1 

1, A 
1, A 
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TABLE 1 



Continued 



Name 


z 


Rmax 


Ls 


L s 


BH mass 


BH mass 


Ref.* 








quasar 1 


quasar 2 


quasar 1 


quasar 2 








(kpc) 


(1O 12 L ) 


(1O 12 L ) 


(10 9 M o ) 


(10 9 M o ) 





2QZJ1217+ 


0006 


1.78 


437.70 


2QZJ1217+ 


0055 


0.90 


315.54 


SDSSJ1226-0112 


0.920 


396.67 


SDSSJ1300-0156 


1.620 


377.02 


2QZJ1328-0157 


2.37 


429.23 


2QZJ1354-0108 


1.99 


464.84 


SDSSJ0846 


f2710 


2.17 


38.42 


SDSSJ0927 


f0707 


1.38 


40.59 


SDSSJ1004 


f4112 


1.73 


31.82 


SDSSJ1007 


f0126 


1.23 


34.85 


SDSSJ1207 


f 1408 


1.77 


33.38 


SDSSJ1431 


f2705 


2.38 


48.21 


SDSSJ1447 


f6327 


1.23 


33.19 


SDSSJ1649 


f 1733 


2.08 


30.17 


SDSSJ1158 


fl235 


0.60 


23.73 


SDSSJ1320+3056 


1.60 


40.16 


SDSSJ1418+2441 


0.57 


29.44 


SDSSJ1426+0719 


1.31 


35.80 


SDSSJ1430+0714 


1.25 


45.14 


SDSSJ1458+5448 


1.91 


43.21 


SDSSJ1606+2900 


0.77 


25.55 


SDSSJ1635+2911 


1.59 


41.68 



2.71 ±0.15 

3.74 ±0.21E-1 
5.46 ±0.30 

1.00±0.06E+1 
7.76 ±0.43 
5.17 ±0.29 
2.15 ±0.12 

9.11 ±0.50E-1 
8.09 ±0.45 

3.80 ±0.21E-1 
1.71 ±0.09 
3.73 ±0.21 

4.02 ±0.22E-1 
6.59 ±0.36 

1.77 ±0.10E-1 
6.83 ±0.38 

2.82 ±0.16E-1 
1.08 ±0.06 
1.69 ±0.09 
1.79 ±0.10 
1.33 ±0.07 
1.23 ±0.07 



4.96 ±0.27 
4.84±0.27E-1 
5.97±0.33E-1 

2.58 ±0.14 
5.08 ±0.28 
2.85 ±0.16 
2.23 ±0.12 

5.70±0.31E-1 
5.53 ±0.31 

6.03±0.33E-1 
1.91 ±0.11 

5.56 ±0.31 

1.57 ±0.09 
4.43 ±0.24 

1.43±0.08 E-l 

1.97 ±0.11 
1.16±0.06E-1 
5.68±0.31E-1 
7.98±0.44E-1 

1.41 ±0.08 
1.21 ±0.07 
5.37 ±0.30 



1.57±3-g? 
2.37±f;|lE-l 
9.75±0.18E-1 
7.84±g||E-l 

A otilS.SO 
4 - JO=t 2.98 
2.94±5-80 

1 - 26±|1 
5.52±i° f f E-l 

A KQ-L8.64 

2.40±fI|E-l 
1 01±n§Q 

2- 14±? : | 
2.53±f-!|E-l 

3.72±*™ 
1.17±g;pE-l 

1.81±3;||E-1 
6.49±f;f V E-l 
q 07+19.00 

i n5+ 2 - 01 
7.91±i 5 J f 



'E-l 



2.82±5-3| 
3.02±|;|E-1 
3.69±X'9?E-1 

|12T- 



1.50± 
2.89± 
1.65± 
1.30± 



k 
k 
\ 

o £.-1+6.94 
O.OO^2.40 ± 

i. 59 4if 

3.73±H2E-1 

1 13p 15 
Q-77 

Q 1 K_I_6.01 

9.32±i 7 . ; rE-l 



c 6.3 

2.53± 



7.88±J 5 . 



3.52±! 



E-2 



7.36± 



5.01 



E-l 
3 E-1 



iE-1 

4.87±||E-1 
8.37±i s 7 9 . 4 E-l 
7.23±Pd| 7 E-l 

Q nK+5.82 



1 
1 

1, A 

1, A 
1 

1 
2 
2 

2, A 
2 

2 
2 
2 
2 
3 
3 
3 
3 
3 
3 
3 
3 



